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We present results of high-resolution thermal-expansion and specific-heat measurements on single 
crystalline a'-NaV20s . We find clear evidence for two almost degenerate phase transitions as- 
sociated with the formation of the dimerized state around 33 K: A sharp first-order transition at 
Ti = (33 ± 0.1) K slightly below the onset of a second-order transition at T% nset w (34 ± 0.1) K. 
The latter is accompanied by pronounced spontaneous strains. Our results are consistent with a 
structural transformation at Ti induced by the incipient spin-Peierls (SP) order parameter above 
T2=Tsp. 



PACS numbers: 75.40.Cx, 75.10.Jm, 65.70.+y 
I. INTRODUCTION 

The discovery of a spin-Peierls (SP) transition in the 
inorganic compound CuGeOs by Hase et al£i has trig- 
gered a large amount of investigations concerning the 
physical characteristics of the SP phase transition. A 
detailed analysis of the magnetic susceptibility showed 
that CuGeOs actually behaves differently from an ideal 
SP system. A significant next-nearest-neighbor coupling 
has to be invokedj-to understand the magnetic properties 
of this compoundH. _ 

Recently, Isobe and Uedatl reported susceptibility results 
on another inorganic system, polycrystalline a'-NaV205, 
which shows a similar rapid decrease in the magnetic sus- 
ceptibility below 34 K. This feature together with the ex- 
istence of a broad maximum in the susceptibility around 
350 K which closely follows the prediction for an S=^ an- 
tiferromagnetic Heisenberg chain led the authors to pro- 
pose a'-NaV20*5 as a possible new SP compound. 
In fact, the isotropic decrease of the magnetic suscep- 
tibility found for single crystalline materialo, the occur- 
rence of a lattice distortion, below the transition tempera- 
ture as observed by X-rayQ, Raman-seatteringp, thermal- 
expansionQ and NMR-measurementsQ as well as the-for- 
mation of a gap in the magnetic-excitation spectrumtl are 
prominent characteristics of a spin-Peierls transition. In 
general, for such a phenomenon to take place, the pres- 
ence of spin- 1/2 chains is considered a necessary con- 
dition. From early structural investigations on polycrys- 
talline material^ one-dimensional (ID) magnetic charac- 
ter was inferredcl, implying two kinds of V chains running 
parallel to the b axis, one of which is formed by magnetic 
(S=l/2) V 4+ , the other one by nonmagnetic V 5+ ions. 
Such a charge-ordered state at T=300K, however, was 
found to be incompatible with recent results of Raman- 
scattering experiments!! Moreover, a refined analysis of 
room temperature X-ray data indicate the existence of 



only one distinct V siteE£l, posing the question whether a 
ID magnetic structure is actually realized in this system. 
So far, most of the work on a'-NaV20s was aiming 
at an investigation of the magnetic and structural as- 
pects whereas little is known about the thermodynamic 
properties. In this paper we will present results of 
thermal-expansion and specific-heat measurements on an 
a'-NaV205 single crystal focussing on the phase transi- 
tion into the nonmagnetic low-T state. 



II. EXPERIMENTAL 

Single crystals of a'-NaV205 were grown by slowly 
cooling a mixture of NaV03 and VO2 with a molar ra- 
tio of 10:1 in an evacuated quartz tube. This assures 
that the melt is always Na rich - a condition necessary to 
achieve a Na content x=l_thc highest possible Na con- 
centration in a'-Na a; V2 05li2l. This is especially important 
as the magnetic properties of a'-Naa;V205 are strongly 
x dependent: already a slight Na deficiency removes the 
SP transition completelya The sample was cooled down 
from 800° C with a cooling rate of 1 K/h. The single crys- 
tals grown by this method could be easily separated from 
the NaVOVflux by boiling the obtained sample in water 
since NaV03 only is soluble in hot water. The suscepti- 
bility of the single crystal (with outer dimensions a x b x 
c = 1.17 x 5.25 x 0.29mm 3 and mass m=5.41mg) selected 
for the present investigations shows a sharp transition at 
Tgp ~ 33 K, accompanied by a rapid and isotropic de- 
crease in x(T) towards lower temperatures (not shown). 
The small upturn in x(T) for T^0 is ascribed to the 
presence of impurities (or defects) . Assuming their effec- 
tive spin to be S=l/2, their concentration is estimated 
to be <0.05%, reflecting the high quality of the crystal. 

The thermal expansion was measured by means of 
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an ultra-high-resolution capacitance dilatometeitia with 
a maximum sensitivity of A1/1=1CP 10 . Length changes 
A1(T)=1(T)-1(4.2K) were detected employing a heating 
rate 4£ < 2.5 K/h. The coefficient of thermal expan- 
sion a(T) = j^-^ was approximated by the differen- 
tial quotient a(T) « p^^g^ with T=(T 1 +T 2 )/2. 
Measurements were performed along the three princi- 
pal axes. The specific heat was measured using an AC- 
calorimeterca. A resistance heater is fed by an AC-power 
with frequency f/2 causing the temperature of the sam- 
ple to oscillate with /. Within an appropriate range of 
frequencies, the amplitude of the temperature oscillation, 
Tac, is inversely proportional to the heat capacity, C(T), 
of the sample. Tac was determined using a lock-in am- 
plifier. The calorimeter was calibrated by measuring a 
sample of pure copper. 

III. THERMAL EXPANSION 

Figure [l] gives a survey of the linear thermal-expansion 
coefficient for T<200K measured along the three princi- 
pal axes. The data reveal a pronounced anisotropy con- 
sistent with the results of an X-ray studycl. While the 
c-axis expansivity is extraordinarily large with a broad 
maximum around 175 K, the b- ('in-chain direction') and 
a-axis expansion coefficients are very small for temper- 
atures 40K<T<250K. A blow-up of the data (cf. inset 
fig.|f|) shows that a a and a& are negative up to about 
150 K yielding a broad minimum around 85 K and 75 K, 
respectively. The low-temperature a(T) behavior is dom- 
inated by sharp phase-transition anomalies around the 
spin-Peierls transition temperature T$p «33K. Figures 
U a-c give details of a(T) near Tsp along the three prin- 
cipal axes (note the different scales). Along the a- and 
c-axis we find large A-type phase-transition anomalies of 
comparable magnitude but opposite signs. A striking 
feature of both data sets is the steepening of the slope of 
a(T) on the low-T side of the transition, i. e. slightly be- 
low the temperature of the a(T) extremum. As becomes 
more evident from the b-axis expansivity data (fig. |^c) we 
assign this behavior to the combined effects of two phase 
transitions, which are almost degenerate in temperature 
(see discussion in section M) : An extremely sharp phase- 
transition anomaly (at Ti) located right at the middle 
of a somewhat broadened one (at T 2 ). While the two 
corresponding phase-transition anomalies are positive for 
a a and negative for a c , they have opposite signs for cxb- 
Moreover, a careful inspection of the a(T) data along the 
three axes shows that the sharp peak is actually split up 
into two very sharp structures separated by only 200 mK. 
This is exemplarily shown for the b-axis data in the inset 
of fig. ^. To exclude the possibility that this fine structure 
is caused by inhomogeneities in the Na content, further 
experiments on an independently prepared single crystal 
are in preparation. 

The influence of a magnetic field of B=8 T applied paral- 



lel to the b-axis on a>b is shown in fig. ||. Surprisingly, the 
field does not lead to a homogenous shift of the two tran- 
sitions towards lower temperatures. Rather it acts quite 
differently on both sides of the transition at Tj.: While 
the high-temperature flank of the second-order transition 
(as e. g. , its onset temperature TJf 156 *) shows a compara- 
tively strong reduction of AT 2 onse * = -(0.36 ± O.f ) K, no 
significant field dependence is found for ctb for tempera- 
tures below Ti, i. e. for T<32.2K. Compared to the effect 
on T^™* 6 *, Ti shows a much weaker field dependence of 
only ATi=-(0.13±0.1)K. 

In order to clarify the nature of the two transitions and to 
separate their respective contributions, we plot in figs. ^ 
and H the relative length changes, AL./1, for i = a, b 
and c on the same temperature scale. As becomes ob- 
vious from fig. ^, the ordering phenomenon below 33 K 
is accompanied by the appearance of pronounced spon- 
taneous strains with large and partly compensating ef- 
fects along the c- and a-axis and an only minor length 
change along the b-axis. Besides this continuous lattice 
distortion, a detailed inspection of the data close to the 
transition provides evidence for the first-order charac- 
ter of the sharp structure in a(T) (cf. fig. ||). This be- 
comes most clear for the b-axis data, where upon lower- 
ing the temperature a sudden elongation occurs within a 
rather small temperature interval of only 0.45 K centered 
around Ti = (33.0 ± 0.1) K. This first-order transition 
is superimposed on the continuous change in the slope 
of Al/l from positive (aj > 0) for T<32.7K to negative 
(a>b < 0) for T>33.25K, which corresponds to the some- 
what broadened second-order transition. We stress that 
in subsequent runs both transitions (including the split- 
ting of the transition at Ti) could be reproduced in detail 
without any significant hysteresis at Ti upon cooling and 
warming through the transition temperature. 
By projecting the finite width of the transition at Ti 
on the a- and c-axis data (vertical lines in fig. ||), one 
can clearly identify first-order-like behavior, i. e. sudden 
length changes, in the Al/l data also for the latter two 
axes. 

For the determination of T 2 , the transition temperature 
of the second-order phase transition, we use the b-axis 
data where, due to their opposite signs, both transi- 
tions can be most clearly separated from each other. 
By simply subtracting in fig. |^c the sharp anomaly at 
Ti from the a&(T) data and replacing the remaining 
broadened transition by an idealized sharp one we find 
T 2 = (32.7 ± O.f) K. We note that a non- vanishing con- 
tribution of the transition at Ti to a(T) at lower tem- 
peratures (neglected in the above procedure) might shift 
T 2 to slightly higher values. 

Our observation of pronounced spontaneous strains be- 
low T 2 indicate a substantial coupling between magnetic 
and lattice degrees of freedom for the present system. 
Within a Landau theory this coupling can be treated by 
adding to the Gibbs free energy an interaction term 
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G int = km, n eTQ n K n G A/"), (1) 

where £i, mj n are the coupling constants between the or- 
der parameter Q an d the strain elements q along the 
ith direction. In general, within the constraints given 
by the symmetry of the system, different coupling types 
(n,m >0) may be possible. For a linear-quadratic strain- 
order-pacameter coupling, ^e^Q 2 , as e.g. , realized in 
CuGeOalij, the condition of a stress-free system then 
yields q(T) oc Q 2 (T), i.e. the spontaneous strains vary 
as the square of the order parameter. To extract the et 
(i=a,b,c), associated with the spin-Peierls transition, for 
the present system one has to subtract from the Aij/Z 
data (cf. fig. H|) the non-relevant background contribu- 
tion. To this end we employed a smooth extrapolation 
of the Al/l data from T>T 2 into the low-T phase. Al- 
though, due to the rather high transition temperature 
substantial uncertainties have to be invoked for the back- 
ground contribution, an only minor effect is found as 
for the conclusion drawn below. The resulting sponta- 
neous strains normalized to their respective values at 4 K 
are plotted in fig. |^ as a function of reduced tempera- 
ture T2 ~ T on logarithmic scales. For all three axes we 

find power-law behavior, oc ( T j^ T ) , over almost one 

decade in temperature. Deviations at low temperatures 
are likely to be caused by uncertainties in the background 
contribution. The tendency towards smaller exponents 
visible for temperatures close to the transition is assigned 
to the influence of one-dimensional fluctuations. From 
the slopes of the linear parts in fig.^ the exponent 7 
can be determined. Most surprisingly, the 7-values for 
the a-, b- and c-axis strains of respectively (0.23±0.05), 
(0.7±0.15) and (0.36±0.1) differ significantly from each 
other. In particular 7^, the exponent for the 'in-chain' 
strain is strongly enhanced - the same tendency, although 
much less pronounced was found also for CuGeOslii As 
for the origin of this behavior one may think of (i) an 
anisotropy in the coupling constants £j,m,n which causes 
different exponents in the leading terms of eqn. [j] along 
the various axes or (ii) a finite coupling to the order pa- 
rameter of the first-order transition. Since for the b-axis, 
the a(T) anomalies associated with the respective tran- 
sitions at Ti and T2 have opposite signs, the one at T2 
being comparatively small, a rather strong influence of 
the lattice distortion at Ti on e(, appears plausible. 

A further observation made in the present thermal- 
expansion study is worth mentioning: A comparison of 
the volume-expansion coefficient, /3, as determined from 
the present measurements on single crystalline material, 
(3 SC = a a + at, + a c , with that derived from polycrys- 
talline samples, /3 pc = 3 • a pc , reveals striking discrepan- 
cies (see ref. || for a pc ). Although Tsp is almost identical 
for the two variants, the absolute value as well as the 
temperature variation of /3 is strikingly different both 
at and above Tgp. A similar observation was already 



made for CuGeC/3 (see ref. [Ti] and [l5| for poly- and single- 
crystalline samples, respectively). Up to now we have no 
explanation for this discrepancy. Texture effects in the 
polycrystal used in refJIJ at the origin of this problem 
can be safely excludedll 

IV. SPECIFIC HEAT 

In order to give a full thermodynamical analysis of 
the two coexisting phase transitions, we have performed 
supplementary measurements of the specific heat, C(T), 
on the same single crystal in the temperature range 
1.8K< T < 60 K. Fig. ^a shows the specific heat as 
C(T) vs T over the whole temperature range investi- 
gated. The SP transition manifests itself in the sharp 
maximum around 33 K. An excellent fit of the low- 
temperature data, i.e. T<12K, is provided by adding 
to the phonon term C p h = /3 • T 3 a magnetic con- 
tribution C mag = aexp(— A/T), which corresponds to 
singlet-triplet excitations across an energy gap A. For 
the zero-temperature value we find A = (84 ± 10) K, 
in good agreement with the results from_measurements 
of the magnetic susceptibility (85±15K)ta, the nuclear 
magnetic resonance CQ8 K)0, as well asJJae electron-spin 
resonance (92±20K)El and (85±20K)EL These values 
are slightly smaller than A=114K_dcrived from inelas- 
tic neutron-scattering experimentsEl. Fig. |7fD shows de- 
tails of C(T) around 33 K as C(T)/T vs T. With the 
knowledge of a(T) (fig. f| one might also separate two 
different contributions to C(T): A sharp maximum as- 
sociated with the first-order transition at Ti (width of 
shaded area in fig. is taken from the anomaly in otb(T) 
at Ti, cf. fig.||c) superimposed on a somewhat broad- 
ened second-order transition at T 2 . We stress that as a 
consequence of the AC technique employed here imply- 
ing temperature oscillations of the sample of ±50 mK at 
T=33K, sharp structures become necessarily broadened 
in our C(T) data compared to a(T). 

V. DISCUSSION 

The results presented in the preceding paragraphs sug- 
gest the existence of two closely spaced phase transitions 
associated with the formation of a dimerized state in a' - 
NaV205 around 33 K: A sharp first-order transition at 
Ti = (33.0±0.1) K slightly above a continuous phase tran- 
sition at T 2 =(32.7±0.1)K. 

As an alternative interpretation one might think of the 
possibility that the observed features reflect a single 
phase transition which is of weakly first order. Since 
the b-axis a(T) data around the transition temperature 
reveal two anomalies of opposite signs, this scenario im- 
plies that the sharp but small negative peak in on, must 
be a 'secondary effect' caused - via the three-dimensional 
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elastic couplings - by the much larger effects in the or- 
thogonal directions. As will be discussed below, however, 
this interpretation appears incompatible with the field 
dependence of ctb and is thus discarded. 

The occurrence of spontaneous strains below T2 sug- 
gests that the latter transition is driven by the strain- 
order-parameter coupling resulting in the formation of a 
gap in the spin-excitation spectrum, i. e. T 2 =Tsp. Then 
the questions at hand are: (i) what is the nature of the 
first-order transition and (ii) what is its role in the forma- 
tion of the spin-Peierls ground state? As for the nature of 
the transition at Ti it is tempting to assign this effect to 
a structural transformation. The substantial magneto- 
elastic coupling, as highlighted by the pronounced spon- 
taneous strains below Tgp, then naturally implies a cou- 
pling of both order parameters. 

As a possible scenario we propose that, upon lowering 
the temperature, the incipient spin-Peierls order param- 
eter induces via strain-order-parameter coupling a struc- 
tural transition at Ti. The latter then stabilizes the 
spin-Peierls transition. Note that along the a- and c- 
axis, where the dominant effect on the lattice is ob- 
served upon cooling below 33 K, the length changes as- 
sociated with the transition at Ti, AI\t i: proceed along 
the same direction as e a and e c , respectively (cf. fig.[|). 
The above assignment of the two transitions is also con- 
sistent with the response in ab when a magnetic field 
of 8T is applied along the b-axis (cf. fig. |): While the 
onset of the transition at T2 shows a significant reduc- 
tion of AT 2 °" set =-0.36K, Ti is only slightly reduced in 
this field by ATi=-0.13K. Even more striking is the fact 
that the low-T end of the continuous transition at T2, 
i. e. CKb(T) for T<Ti, does not show any significant field 
dependence at all. These shifts have to be compared 
with the field- induced reduction, ATgp , expected for an 
ideal spin-Peierls transition^-, 
known spin-Peierls systemstBu. 
a Hartree-Fock approximation 



tions would be provided by pressure experiments, where 
uniaxial strain is applied along the b-axis. From the dif- 
ferent signs of the phase-transition anomalies Aab and 
Alb associated with the respective transitions at Ti and 
T2 , a pressure-induced separation of the transitions can 
be expected: Applying the Ehrenfest-relation and the 
Clausius-Clapeyron equation to the second- and first- 
order transition, respectively, we predict a positive pres- 
sure dependence of the transition at T2 ( ( |^ ) > 0) 
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and a negative pressure dependence for the first-order 
transition ((gi)^ < ). 

In conclusion, thermodynamical investigations on a'- 
NaV 2 05 reveal indications for two almost degenerate 
phase transitions around 33 K associated with the for- 
mation of a nonmagnetic ground state. The second- 
order transition at T2=(32.7±0.1) K shows the signatures 
of a SP transition as to the occurence of spontaneous 
strains and the field dependence of its onset tempera- 
ture. In this sense we assign the first-order transition at 
Ti = (33±0.1)K to a structural transformation induced 
by the incipient SP order parameter at T>T2. In general, 
for this phenomenon to take place, the existence of spin- 
1/2 chains is considered a necessary condition. In light 
of the results of recent structural investigations at room 
temperature which are inconsistent with the above re- 
quirements a detailed T-dependcnt determination of the 
magnetic structure appears badly needed. 
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,and actually found for all 
Employing the results of 
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where g&=1.972 is the g-factorEH, /i# the Bohr magne- 
ton, kg the Boltzmann constant and a=0.44Ea, we find 
ATgp =-0.375 K. From the excellent agreement with the 
experimentally observed reduction of T° nset we infer that 
the high-T side of the transition, i. e. at T>Ti, reflects 
the behavior of a common spin-Peierls transition. The 
much weaker field dependence of Ti and the almost field- 
independent behavior below Ti rule out the possibility 
that the sharp feature in ab at Ti is a secondary effect 
caused by the pronounced anomalies in the orthogonal 
directions. In this case a homogenous shift of the two 
features at Ti and T2 in a magnetic field would be ex- 
pected. Rather it indicates the action of another inter- 
fering mechanism. 

A direct proof of our assignment of two phase transi- 
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FIG. 7. Specific heat, C(T), as C vs T (a) and C(T)/T vs 
T (b) . Width of shaded area in (b) indicates the temperature 
range of the first-order transition as determined from a(T) 
measurements. 



FIG. 1. Coefficient of thermal expansion, a, vs T along the 
three orthorhombic crystal axes. Inset: blow up of a vs T for 
the a- and b-axis. 



FIG. 2. Coefficient of thermal expansion, a, vs T for 
28K<T<37K along the various axes. Inset shows details of 
the negative peak in ctb on expanded scales. 



FIG. 3. Coefficient of thermal expansion, at, for B=0 (□) 
and B=8 T (•), applied parallel to the b-axis. Lines are guides 
to the eyes. Data points around the double structure near Ti 
have been omitted for clarity. 



FIG. 4. Relative length changes Al/l vs T along the three 
crystal axes shifted so that the curves meet at 35 K. 



FIG. 5. Relative length change Al/l vs T on expanded 
scales close to the SP transition. Vertical lines indicate the 
temperature range of the first-order transition. 



FIG. 6. Normalized spontaneous strains vs reduced tem- 
perature on logarithmic scales. Solid lines represent fits of 
the form e cx (1 -T/T 2 ) 7 . 
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